The eccentricity in heavy-ion collisions from Color Glass Condensate initial conditions 
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The eccentricity in coordinate-space at midrapidity of the overlap zone in high-energy heavy- 
ion collisions predicted by the fcx-factorization formalism is generically larger than expected from 
scaling with the number of participants. We provide a simple qualitative explanation of the effect 
which shows that it is not caused predominantly by edge effects. We also show that it is quite 
insensitive to "details" of the unintegrated gluon distribution functions such as the presence of 
leading-twist shadowing and of an extended geometric scaling window. The larger eccentricity 
increases the azimuthal asymmetry of high transverse momentum particles. Finally, we point out 
that the longitudinal structure of the Color Glass Condensate initial condition for hydrodynamics 
away from midrapidity is non-trivial but requires understanding of large- a: effects. 

PACS numbers: 12.38.Mh,24.85.+p,25.75.Ld,25.75.-q 

I. INTRODUCTION 

Elliptic flow [1] - the azimuthal momentum anisotropy 
of produced particles with respect to the reaction plane, 
i>2 = (cos(2</>)) - is one of the key observables for collec- 
tive behavior in heavy ion collisions. The large V2 ob- 
served experimentally in semi-central Au+Au collisions 
at the Relativistic Heavy-Ion Collider (RHIC) [2] is con- 
sistent with non-viscous hydrodynamic expansion of a 
quark gluon plasma (QGP) droplet [3]. On the other 
hand, a hadron "plasma" with small transverse pressure 
as initial condition underestimates 1*2 at midrapidity [4]. 
These observations have been interpreted as an indica- 
tion for the formation of a QGP with very small viscosity 
shortly after the instant of collision. 

Specific solutions of partial differential equations such 
as hydrodynamics in general depend on the assumed 
boundary conditions. Here, we focus on the initial condi- 
tions for the subsequent evolution of the matter created 
in non-central Au+Au collisions at RHIC. Specifically, 
we consider the "Color Glass Condensate" (CGC) initial 
conditions within a fc^-factorization approach [5]. 

The elliptic flow vi in the final state is proportional to 
the initial spatial eccentricity e in the plane r± = (r x , r y ) 
perpendicular to the beam axis, 



in the cool hadronic phase are taken into account (by us- 
ing a hadronic cascade [7] rather than hydrodynamics). 
However, ref. [6] argued that the different eccentricities 
obtained from the CGC and Glauber model, respectively, 
arise from mainly the low-density edges of the overlap re- 
gion where the CGC approach is less trustworthy. 

In this paper we show that more sophisticated models 
for the nuclear unintegrated gluon distribution (j)(x,k±) 
than those employed in the original KLN approach [5] 
and in [6] also predict much larger eccentricity than the 
Glauber approach. In particular, we shall calculate e us- 
ing gluon distributions which account not only for satura- 
tion but also for extended geometric scaling (anomalous 
dimension) and leading-twist shadowing, both of which 
are important for transverse momenta above the satu- 
ration momentum Q s . We also provide a simple quali- 
tative explanation of the large eccentricity in the CGC 
approach which motivates that it does not originate from 
pure "edge effects". We therefore expect that both the 
large eccentricity of the transverse overlap region as well 
as the in-plane "twist" of produced gluons [8] are rather 
generic predictions of the k± -factorization approach. Fi- 
nally, using a simple model for energy loss induced by the 
dense medium, we show how high-pj_ partons could probe 
the initial distribution of produced gluons in coordinate 
space. 



(1) 



The average is taken with respect to the initial energy 
density distribution. The reaction plane is spanned by 
the beam (z) and impact parameter (x) directions. 

Ref. [6] observed that the CGC initial conditions cor- 
respond to much larger eccentricity than the commonly 
assumed scaling of the initial energy density with the 
local number of participants. In fact, hydrodynamical 
simulations with such initial conditions were shown to 
overestimate the measured vi even if dissipative effects 



II. GLUON PRODUCTION 

In the fc^-factorization approach the number distribu- 
tion of produced gluons is given by 



d?r±dy 



an c r p ± d 2 P± 

X 4>b(x 2 , (p±-fc±) 2 



d 2 k± a s (j>A{xi,Px) 
(2) 
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with N c = 3 the number of colors. Here. p± and y de- 
note the transverse momentum and the rapidity of the 
produced gluons, respectively. The light-cone momen- 
tum fractions of the colliding gluon ladders are then given 
by £1,2 = P_l exp(±y)/ v /ijv r /v, where ^/Jnn denotes the 
center of mass energy. We set p™ ax such that the mini- 
mal saturation scale Q" un (xi) in the above integration is 
h-QCD = 0.2 GcV. We note that the eccentricity at midra- 
pidity is rather insensitive to the choice of p™ ax [9] . The 
KLN approach [5] employs the following unintegrated 
gluon distribution function at x <C 1: 



Due to < 7 < 1, §{k\ > Q 2 ) from (6) is harder 
than the one from the original KLN model. The con- 
tribution from the extended geometric scaling regime to 
the p^-integrated multiplicity is small when Q s is large. 
Since the eccentricity e involves an integration over a 
wide range of saturation momenta, however, it appears 
less obvious whether or not the geometric scaling window 
above Q s plays a role. Here, we employ the following pa- 
rameterization of 7 for k± > Q s [13]: 



7(xi,fe|) 



Is + (I -Is) 



log(fci/<&) 



1 



XYi + 1.2^ + ]og(ki/Qi i ) 



a s (Q1) max(Q2,fc2) 



(3) 



Q s denotes the saturation momentum at the given mo- 
mentum fraction x and transverse position rj_ . The over- 
all normalization of 4> is determined by the multiplicity at 
midrapidity for the most central collisions. In the KLN 
approach, the saturation scale for a nucleus in A+B col- 
lisions is then parameterized as 



Q^, B (x,r ± ) = 2GeV 2 




(4) 



tpart(^j_) is the transverse density of participants. 



where n p 

It is calculated from the thickness functions Ta and Tb'- 



x (l-(l-a NN T B (r ± -b/2)/B) B ) .(5) 

b = (6, 0) is the impact parameter vector in the trans- 
verse plane. The nucleon-nucleon inelastic cross section 
at ^/sjvjv = 200 GeV is <jnn — 42 mb. 



(7) 

where 7 S = 0.627 is the anomalous dimension corre- 
sponding to leading-order (LO) BFKL evolution with 
saturation boundary conditions in the saddle point ap- 
proximation [12]. The rapidities Yi of the fusing gluon 
ladders are related to their momentum fractions x\ by 
Yi = ]ogQ./xi). 

The extended geometric scaling window corresponds 
to the region of k± where 7 evolves from 7 S to ~ 1. For 
asymptotically small x, one can neglect the term ~ s/Y 
in (7) and the upper end of the extended geometric scal- 
ing window is then given by Q gs {x) ~ Q 2 (x)/Qo [12]; for 
the central region of RHIC collisions, the approach to the 
DGLAP regime is mainly driven by the subasymptotic 
~ y/Y term. The above form of 7 reproduces the evo- 
lution of inclusive ^-distributions from minimum-bias 
d + Au collisions with rapidity rather well [13]. 

Finally, we shall also compute e using an unintegrated 
gluon distribution obtained from an Ansatz for the dipole 
forward scattering amplitude N{x,u±) [14, 15]: 



The form Qj(x) 



.-A 



with A « 0.2 — 0.3 is motivated 



0, k 2 ± ;r ± ) 



1 



d 2 u 



— ih\ -u 1 rj2 



ViN(x,u x ), (8) 



where u± is the transverse size of the dipole, conjugate 
to k±. A viable model for the unitarized dipole profile 
is [13, 14] 



by deep inelastic scattering (DIS) data from the Hadron 
Electron Ring Accelerator (HERA) for x < 0.01 [10]. 
The parameterization (4) leads to an average saturation 
momentum at midrapidity of about 1.5 GeV for central 
Au+Au collisions at ^/snn = 200 GcV, which is in line 
with standard estimates [5]. 

The KLN model (3) exhibits a transition from satu- 
ration, ~ const., directly to the DGLAP regime wne re the anomalous dimension 7(2;, uV) is of the same 



N(x, u±) = 1 — exp 



1 



2 n 2\ 7 (x,ul) 



(9) 



0(/fj_) ~ a * fc-L = Qs', it docs not incorporate 

leading-twist shadowing since <f>(k± > Q s ) is proportional 
to the density of participants at a given transverse coor- 
dinate. On the other hand, it is well known that DIS 
data from HERA exhibit approximate "geometric scal- 
ing" even for Q 2 > Q 2 [11]. This indicates that at small 
x the anomalous dimension 7(fc_i_) of the gluon distribu- 
tion function evolves smoothly towards its DGLAP limit 
of 7 = 1 — 0(a s ) ~ 1 as k±/Q s increases. This is a con- 
sequence of quantum evolution of </> with log 1/x [12]. 
The unintegrated gluon distribution can then be written 
as 



(x,k 2 j_;rj_) 



s(Q 2 ) 



min 




7(a,fci)N 



(6) 



form as (7) but with the replacement kj_ — > l/uj_. This 
dipole profile describes the transitions from saturation 
(N(u±) ~ 1) to the leading- twist extended geomet- 
ric scaling (N(u±) ~ (u^Q 2 ) 1 ) and DGLAP regimes 
(N(u±) ~ u^Q 2 ). Leading-twist shadowing is present 
in the extended geometric scaling window since with 
Q 2 ~ A 1 / 3 , N ~ A 7 / 3 exhibits a weaker dependence 
on the nuclear thickness. 



III. RESULTS 

In this section, we show the results of a numerical inte- 
gration of eq. (2), using the above-mentioned models for 
the unintegrated gluon distribution cj>(x,k±). We shall 



3 



5! 4 
c 



; \^\ KLN,A,=0.2~0.3 


! (KLN) y .X=0.2~0.3 




200 ^jjf^ 








130 




, i , , , , i , , , 



100 



200 

Npart 



300 



400 



CM 



KLN,X=0.2 ~ 0.3 
(KLN) Y , A,=0.2 
DHJ Dipole 




200 



100 



200 

Npart 



300 



400 



FIG. 1: (Color online) Centrality dependence of the multiplic- 
ity in Au+Au collisions as compared to PHOBOS data [17] 
for y/SNN = 130 and 200 GeV. Running (fixed) coupling is 
used in the upper (lower) panel. A denotes the growth rate 
of the saturation momentum with log 1/x. The curves corre- 
spond to various parameterizations of the unintegrated gluon 
distribution function, see text for detailed explanations. 



employ both fixed (a s = 0.5) and LO running coupling 
constants, a s (Q 2 ) = 4n/9\og(Q 2 /A 2 QCD ) [16]. 

We first check the centrality dependence of the multi- 
plicity. In Fig. 1 we compare the calculated number of 
gluons to PHOBOS data [17] for charged hadrons, as- 
suming that the hadron yield is proportional to the yield 
of produced gluons. "(KLN)" denotes the original KLN 
parameterization of the unintegrated gluon distribution 
function (3), "(KLN) 7 " its extension from eq. (6) to in- 
clude the smooth transition of the anomalous dimension 
and leading-twist shadowing, and "DHJ Dipole" the uni- 
tarized dipole profile from eqs. (8,9). We find that all 
models are consistent with the data within error bars. 
As a side-remark, we note that the explicit numerical in- 
tegration of eq. (2) over r± improves the KLN results at 
midrapidity, especially in the high iV part region, as com- 
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FIG. 2: (Color online) Initial spatial eccentricity e at 
midrapidity as a function of the number of participants for 
200 A GeV Au+Au collisions from various models. For com- 
parison, we also show initial conditions where the initial par- 
ton density at midrapidity scales with the transverse density 
of wounded nucleons (full line) and of binary collisions (dotted 
line) [19]. 



pared to the Q 2 (r±) — > (Q 2 S ) approximation employed 
in Refs. [5]. We refer to Ref. [9] for a nice explanation 
of the positive slope of the curve; also, to ref. [18] for a 
computation of the centrality dependence of the multi- 
plicity within the classical McLerran-Vcnugopalan model 
and its relation to the KLN approach. 

Next, in Fig. 2, we compare the eccentricity for the 
above-mentioned unintegrated gluon distribution func- 
tions. Regardless of running or fixed coupling, all sat- 
uration models yield almost the same £ as a function of 
centrality while that predicted by Glauber-type models 
is much lower. The latter is commonly used as initial 
condition for hydrodynamical simulations [19], assuming 
that the energy density in the transverse plane scales with 
either the number of participants or the number of col- 
lisions (dominance of hard processes), respectively. The 
fact that all of the above models for <j> predict about the 
same e already indicates that the large anisotropy in co- 
ordinate space might be generic to the fcj_-factorization 
formula (2). 

The large predicted eccentricity could be understood 
qualitatively in the following way. In the fcj_-factorization 
approach the number of produced gluons scales approxi- 
mately with the smaller of the saturation scales of the two 
nuclei [5] . This is valid not only as a function of rapidity 
but also as a function of the transverse coordinate [20] : 



dN 



d 2 r±dy 



min(Q 2 s Jy,r ± ),Ql 2 (y,r ± )). 



(10) 



To understand this, approximate the local unintegrated 
gluon distribution by 4>A(k±) ~ 0(<3s.mm — k±). Seal- 



FIG. 3: Along the line l x , the density of gluons in the CGC 
case falls off more rapidly than in Glauber type models. Along 
l y , the collision is symmetric and the CGC gluon density be- 
haves similar to the Glauber model. See text for explanation. 



ing with Q^min follows also from analytical solutions of 
the classical Yang-Mills equations for asymmetric colli- 
sions [21]. 

Fig. 3 depicts two paths away from the center of the 
overlap region. Along l x the number density scales, to 
logarithmic accuracy, as 

PCGc(r x ,0) ~ Q 2 s l (r x ,0) ~ n part ,i(r x , 0) . (11) 

On the other hand, in the opposite direction along —l x 
we have 

PCGc(r*,0) ~ Q* 2 (r x ,0) ~ n part , 2 (r x , 0) . (12) 

In the wounded nucleon model, the density scales with 
the total number of participants from both nuclei, 

PGlauber^x, 0) ~ U 

part, as? 0) ^part,2 (j'x j 

0) , (13) 

which of course drops less rapidly with r x since the drop 
of the density towards the edge of the more dilute nucleus 
is at least partly compensated by the increasing number 
of participants from the other nucleus. 

Along the line H y , in turn, both saturation scales arc 
equal and so 

/0cgc(O, r y ) ~ (n P art,i(0,r y ) + n par t,2(0, r y ))/2 

~ pGlauber(0, r y ) . (14) 

Since the number density along l x drops more rapidly for 
the CGC initial conditions than for the Glauber model, 
while both behave the same along l y , the resulting eccen- 
tricity is higher for the former. Thus, the large eccentric- 
ity is a feature inherent to the fc^-factorization formula. 

Contrary to simple geometrical overlap models, the 
CGC also predicts a non-trivial longitudinal structure 
of the initial condition. Even under the simplest as- 
sumption of local quantum evolution, which implies 
that the r_L-dcpendence of 4>a,b arises exclusively via 

Ql A B( x ' r ^) ~ n part( r -i-) a t an Y rapidity, moments of 
the density of produced gluons in coordinate space do 
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FIG. 4: (Color online) Moments ((r x - (r x )) n ) (in fin") of the 
coordinate-space distribution of produced gluons as a func- 
tion of rapidity. Open (closed) symbols correspond to the 
Ansatz (3) with (without) the additional factor of (1 — x) 4 . 
The n = 1, 3 moments have been scaled up by a factor of 10. 



depend on y ! For example, the "center of gravity" (r x ) 
does not vanish at y ^ 0, as already discussed in [8]. 

For gluon production at non-zero rapidity the momen- 
tum fraction x of one of the fusing gluon ladders grows 
like ~ exp(y). In the KLN approach, large- a; effects are 
incorporated into the unintcgratcd gluon distribution by 
the substitution <f>(x, k\) — ► <f>(x, k\){\ — x) A . This is, of 
course, a rather qualitative model of large-a; effects which 
does not treat their Q 2 dependence etc. 

Fig. 4 shows the evolution of various moments of the 
energy density of produced gluons with rapidity. The 
longitudinal structure is clearly non-trivial since odd mo- 
ments of r x no longer vanish at y ^ 0. Hydrodynam- 
ical solutions will exhibit a "left-right" asymmetry (i.e. 
vi ^ 0) due to the fact that ((r x ~(r x )) 3 ) ^ [22]. Hence, 
as compared to the geometry at y — illustrated in 
Fig. 3, not only does the center of the "overlap" shift with 
rapidity but the entire shape changes. Fig. 4 also shows 
that the (1 — a;) 4 factor, which models large-a; effects, 
strongly affects the third moment of r x at y> 2. (Other 
"details" such as p™ ax from eq. (2) are also much more 
important at large y, which is another reason why there 
the numbers from Fig. 4 can not be trusted quantita- 
tively.) This result indicates that our present limited un- 
derstanding of "large-a;" effects within the CGC formal- 
ism inhibits more quantitative calculations of the initial 
conditions for hydrodynamic models over a broad range 
of y [23] . In order to compute the observable momentum- 
space distributions of particles from these initial condi- 
tions one needs to solve for the subsequent thermalization 
and three-dimensional hydrodynamical expansion [24] of 
produced gluons, which is beyond the scope of the present 
note. 

However, high-px particles may probe the initial den- 
sity distribution via final-state interactions such as en- 
ergy loss. We employ a simple model for radiative energy 
loss of fast partons to illustrate the qualitative effect. 
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In a dense QCD medium the induced radiative en- 
ergy loss reduces the initial transverse momentum p^_ of 
a produced hard parton prior to its fragmentation into 
hadrons. The "quenching factor" (1 — px) depends in 
general on the initial parton rapidity y, its initial trans- 
verse momentum p^_ as well as on the local density [25] . 
Here, we consider only mid-rapidity (y = 0) partons and 
further simplify the calculation by neglecting the mild 
p± dependence of the fractional energy loss. We use the 
model proposed in [26] and further developed in [27]. p 
is a parameter set to reproduce the observed average jet 
quenching at b = 0. For eikonal (straight-line) jet trajec- 
tories, the opacity x of the medium is a function of the 
initial production point r j_ of the jet and of its azimuthal 
direction of propagation v±_ = (cos <p, sin tp), 

/■oo I — I 
dl — — - p{r ± + lv±;b) . (15) 

where p(r±,b) denotes the local bulk density. Here, we 
assume a formation time of the medium of Iq = 0.2 fm. 
Larger formation times tend to increase the final az- 
imuthal asymmetry [28]. Eq. (15) accounts for non- 
dissipative one-dimensional expansion of the medium 
along the beam axis. 

For simplicity, we approximate the initial p±- 
distribution of hard partons with p° ± > Q gs by a power 

law dcr/d 2 p L dy ~ (p < ±)~ n ^-\ where n(p ± ) is a slowly 
varying function which we calculate in the DGLAP 
leading-twist approximation. This is justified since, by 
definition, the anomalous dimension 7 — > 7dglap m the 
high-px regime. The final transverse momentum of a 
given parton then is p± = p\(l — We define the 

quantity Raa(p.l, <P',b) as the ratio of the initial to the 
final transverse momentum distributions, omitting frag- 
mentation into hadrons: 

R AA ( P± ,p;b) = (f g (p±)(l - fixT sip±) - 2 

+ (l-/ 9 ( p± ))(l-i MX )«^)-2\ > j(16) 

where f g (p±) is the fraction of gluons and n q / g (j>±) arc 
the quark and gluon power laws, respectively [29]. The 
color factor C F /C A = (N 2 - l)/(2N 2 ) = 4/9 accounts 
for the weaker coupling of quark jets to the medium. The 
average in (16) is performed over the initial production 
points r± of the jets. In the DGLAP regime, these are 
distributed according to the local number of binary col- 
lisions, ~ T A (r ± + 6/2) T B {r ± - 6/2). 

Fig. 5 displays the angular average R AA (p±; b) of (16) 
as well as its second harmonic, denoted as l>2(p_U b). The 
former is rather insensitive to the shape of the "overlap" 
region as it involves an integral over the azimuth tp. The 
asymmetry 1)2 for semi-central collisions is slightly larger 
for CGC versus Glauber initial conditions. Nevertheless, 
for the assumed non-dissipative expansion, which leads 
to a rather rapid cooling of the medium, the difference is 
moderate. 
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FIG. 5: Top: suppression factor 7?aa(p±;6) at p± ~ 10 GeV 
as a function of centrality for CGC/KLN and ~ iVpart bulk 
density distributions, respectively. PHENIX data [30] corre- 
spond to hadron transverse momenta p± > 4.5 GeV. Bottom: 
second harmonic coefficient of Raa(p±, tp; b) / Raa(p±; b). 



IV. CONCLUSIONS 

We have shown that the eccentricity predicted by CGC 
models is larger than that from Glauber models. This 
property is due to the different scaling behavior of the 
density in the transverse plane of the collision. Namely, 
k± -factorization predicts that dN/d 2 r± is proportional 
to the smaller rather than to the average of the two local 
saturation scales. We also find that the centrality depen- 
dence of the multiplicity and of the eccentricity depend 
only weakly on "details" of the unintegratcd gluon distri- 
bution function, such as leading-twist shadowing and the 
anomalous dimension. The centrality dependence of the 
saturation momentum can be tested further by high-p^ 
particle production in p+A collisions [31]. 

These observations could be very important for our 
understanding of the collective expansion of the high- 
density matter created at RHIC and, in the near fu- 
ture, at the LHC. Larger initial eccentricity should lead 
to larger V2 in the final state unless the build-up of az- 
imuthally asymmetric collective flow is damped by dis- 
sipative effects. Ref. [6] found that dissipation in the 
hadronic stage might not be sufficient to explain experi- 
mental data. Dissipative corrections to relativistic ideal 
fluid expansion of a QGP have recently been discussed in 
refs. [32] . Previous parton cascade simulations with only 
clastic two-body scattering concluded that a huge trans- 
port opacity (large parton cross sections and/or large 
parton densities) is needed to account for the measured 
V2 [33]. It will be most interesting to reconsider those 
results using the CGC initial conditions suggested here 
together with parton multiplication processes [34]. 

Finally, we point out that a better understanding 
of large-x effects within the CGC is important to pro- 
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vide quantitative results on the structure in coordinate 
space of the initial condition for hydrodynamic models of 
heavy- ion collisions at RHIC and LHC. 



Note added: This note has been added in compliance 
with a request of the referee of this paper to clarify our 
use of the participant density (5) in the definition (4) 
of the saturation momentum, and in the unintegrated 
gluon distribution function (3). This may raise concern 
as the unintegrated gluon distribution functions 4>a and 
4>b then appear to be non-universal quantities which de- 
pend on properties of the collision partner. One might 
prefer to take Q 2 AB to be proportional to the thickness 
function T AB of the corresponding nucleus rather than 
to its participant density. 

We have found that doing so leads to a significant 
overestimate of the multiplicity in peripheral collisions 
and hence to a less good description of the data shown 
in Fig. 1 than for the original KLN approach [5] with 
Q 2 ~ Mpart- This failure can be corrected by taking 
into account explicitly the probabilities pa.b < 1 of 
actually encountering at least one nucleon in nucleus 
A or B, respectively, at a given transverse coordinate. 
This redefines the unintegrated gluon distributions as 
&a.b = Pa.b <t>A,B- The r^-dependence of Q 2 can then 
be defined as Ta,bIpa.b- : which prevents it from drop- 
ping below the saturation scale for a single nucleon at 
the periphery of the nuclei. Factorization is now mani- 



festly preserved and we also find that the ccntrality de- 
pendence of the multiplicity is very similar to that from 
the original KLN approach, shown in Fig. 1. 

With respect to the multiplicity N = dN/d 2 r±dy, 
however, the approach sketched in the previous para- 
graph is essentially identical to that of rcf. [5]. To re- 
alize this, note that the multiplicity is, to a good ap- 
proximation, a homogeneous function of first degree: 
N(kQ 2 a ,kQ 2 sB ) ~ K m N{Q 2 A ,Q 2 sB ), with m = 1. This 
can be easily understood from the approximate form N ~ 
mm(Q 2 A , Q 2 B ) | log Q 2 a /Q 2 b\i f° r example. Since now 
PaPb Q 2 s b i s proportional to paTb = ^p a rt> we recover 

PAPB N(T A /PA, Tb/pb) ^ ^V(npart> n part)- 

The transverse energy dE±/d 2 r±dy ~ 
mm(Q 2 sA ,Q 2 sB ) max(Q sA , Q s b) | log Q 2 aA /Q 2 sB \, 
which we use as a weight in eq. (1), is homogeneous 
to degree m = 3/2 rather than 1. Consequently, we 
find that the approach sketched in this note predicts an 
eccentricity roughly half-way inbetween the Glauber and 
the KLN approaches for the most peripheral collisions 
but is close to the latter for semi-central collisions where 
T Aj b and n p ^ rt differ very little. 
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FIG. 6: Additional Figure: Eccentricity as a function of 
A^part for the KLN-model of the unintegrated gluon distri- 
bution function. The cut-off in the integrations over k± 
and p± in eq. (2) is varied from our "standard" choice of 
Qs{xi) > Aqcd to Qs(xi) > 3Aqcd. Note that a larger 
cut-off effectively chops off part of the low-density rim of the 
"overlap" region. The eccentricity then increases slightly. 
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FIG. 7: Additional Figure: The mid-rapidity energy den- 
sity profile for a Au+Au collision at b = 8 fm as obtained 
from eq. (2) with the KLN unintegrated gluon distribution or 
from a Glauber model, respectively. 



